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Abstract. We introduce the concept of initial-phase spectroscopy as a control of 
the dynamics of entangled states encoded into a two-atom system interacting with a 
broadband squeezed vacuum field. We illustrate our considerations by examining the 
transient spectrum of the field emitted by two systems, the small sample (Dicke) and 
the spatially extended (non-Dicke) models. It is found that the shape of the spectral 
components depends crucially on the relative phase between the initial entangled state 
and the squeezed field. We follow the temporal evolution of the spectrum and show that 
depending on the relative phase a hole burning can occur in one of the two spectral 
lines. We compare the transient behavior of the spectrum with the time evolution 
of the initial entanglement and find that the hole burning can be interpreted as a 
manifestation of the phenomenon of entanglement sudden death. In addition, we find 
that in the case of the non-Dicke model, the collective damping rate may act like an 
artificial tweezer that rotates the phase of the squeezed field. 
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1. Introduction 

Recently there have been considerable interest in practical applications of squeezed 
light. The interest stems from the recognition and experimental demonstrations that 
squeezed light could be used as a resource for entangling atomic ensembles [HE], noise- 
free information transfer and processing of entanglement [Sill]. Testimony to this ability 
is provided by experimental observations of such effects as quantum teleportation [5l [6] , 
quantum erasing [Tj E], quantum dense coding P [TOl [H], creation of entangled cluster 
states [121 IISI Ull US] , and entanglement storage in atomic ensembles [161 113 [H] • 

The early studies of squeezed light and its applications showed that atomic 
dynamics and fluorescence can be significantly modified in the presence a squeezed 
vacuum field [19]. The reason is in an unequal phase- dependent redistribution of 
the fluctuations (noise) between the two quadrature components of the field. It 
has been shown that this modification can be a source of many novel effects and 
techniques in quantum optics [201 EB [22]. In particular, the interaction of atoms 
with a squeezed vacuum field can lead to the inhibition of atomic phase decay [23] . 
subnatural linewidths [241 [25] . hole burning and dispersive profiles in the fluorescence 
and absorption spectra [261 EZ], a linear dependence on the intensity of the two- 
photon absorption rate [281 EH EOl |3T1 [32] , transparency and gain without population 
inversion [331 IS]- The fluorescence and absorption spectra were found to exhibit a 
striking dependence on the relative phase of the squeezed vacuum and the coherent 
driving field [351 [SEl [37]. Depending on the phase, the spectral lines can exhibit a 
subnatural or supernatural linewidth. 

Although the phase dependent effects are usually predicted for systems that interact 
with a squeezed vacuum and are simultaneously driven by a coherent field, there can 
be interesting phase- dependent effects present even in the absence of the coherent 
field. This situation might be encountered when a system is initially prepared in a 
superposition state of its energy eigenstates. To put it another way, phase dependent 
effects are expected to vanish whenever the initial state of the system is an eigenstate of 
the energy, because then the emitted photon has no phase. Entangled states of a multi- 
atom system could provide a reference phase for the phase of the squeezed vacuum field. 
With the recent great interest in the dynamics of entanglement in atomic systems, it 
is worth to study phase dependent changes in the time evolution of entangled states of 
two atoms coupled to a squeezed vacuum field. Needless to say, the relative phase could 
be used as a control of the evolution of the entangled system. 

In this paper we consider what may be termed "initial-phase spectroscopy" to 
search for a directly measurable signature of the abrupt loss of entanglement — the 
phenomenon of sudden death of entanglement [38ll39llini[lll[l2lll3l[lll[l5lll6]. We 
study the transient spectrum of the field emitted from two two-level atoms initially 
prepared in a pure (phase dependent) entangled state and spontaneously decaying 
in a squeezed vacuum field. The squeezed field is treated as a correlated phase 
dependent noise environment. It is well known that an entangled state has a phase 
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and nonvanishing correlations and therefore it is interesting to know how these initial 
correlations affect the spectrum when the system is subjected to a correlated phase 
dependent noise. We illustrate our considerations by examining two different models of 
the two-atom system, in which phase dependent effects arise dramatically. These are the 
small sample (Dicke) model [17] and the spatially extended (non-Dicke) model |1HI Il9] . 
In both cases, the proposed initial-phase spectroscopy is found to be able to detect 
the phenomenon of entanglement sudden death by observing changes in the transient 
spectrum. These changes would also reveal other novel features that could provide 
a measure of the relative phase between two oscillating atomic dipoles. The transient 
spectra are produced in the following manner: in the time period before t = the atoms 
are prepared in an entangled state. At some time {t = 0) following the preparation of 
the atoms in the entangled state, a broadband squeezed vacuum field is applied to the 
atoms. Then, we trace the time evolution of the fluorescence spectrum from t = to the 
steady-state distribution by monitoring the spectrum as a function of time for different 
relative phases of the initial entangled state and the applied squeezed vacuum field. 

We interpret the results in terms of populations and two-photon coherence between 
collective states of the system and derive an analytical formula of the spectrum in 
the limit of the well separated spectral lines. It is well known that the intensity of a 
spectral line corresponding to transition between two collective states is proportional to 
the population of the state from which the transition originates. We find the intensity of 
one of the two spectral lines depends not only on the population of the upper state but 
also on the two-photon coherence. This dependence introduces the possibility to obtain 
the phase sensitive magnitude of this spectral line. In the course of the calculations, we 
observe that the sensitivity to the relative phase has an additional feature characteristic 
of entangled properties of the system. We find that the hole burning in the spectral 
line occurs at time when the atoms become disentangled. We are therefore led to the 
conclusion that observations of the hole burning in the spectral line provide a measurable 
signature of the phenomenon of entanglement sudden death. 

The paper is organized as follows. We start by introducing the system and the 
definition of the transient spectrum. In section [3l we derive the equations of motion 
for the density matrix elements which we then solve for two different cases: the small 
sample (Dicke) model corresponding to atoms confined to a region of dimensions much 
smaller than the resonant wavelength, in section IH and for a spatially extended (non- 
Dicek) model where the atoms are confined to a region of dimensions comparable and 
larger than the resonant wavelength, in section [51 Finally, in section [6], we summarize 
our results. 

2. Transient fluorescence spectrum 

We propose to study phase dependent effects of an entangled system in terms of the 
transient spectrum of the fluorescence field emitted by a two-atom system interacting 
with a broadband multi-mode squeezed vacuum field. Different definitions of the 
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transient spectrum have been proposed, with their weaknesses and strengths summarized 
in Ref . [50] . As we shall be interested in signatures of the phenomenon of entanglement 
sudden death, we assume that the system starts at t = from a well defined initial 
entangled state and use the time dependent physical spectrum [511 

ft rt 



= 21^ Vr,, / dti / dt2 
i,3=i -^0 Jo 



where t is the elapsed time after the system was prepared in the entangled state, 2Vd is 
the bandwidth of the detector which has its transmission peak at frequency cj, S^ and 
S~ are the dipole raising and lowering operators, respectively, of the ith. atom, Va = T 
are the spontaneous emission decay rates of the atoms, equal to the Einstein A coefficient 
for spontaneous emission, and Tij = Tji, {i ^ j) is the collective damping which results 
from an exchange of photons between the atoms. The collective damping depends on 
the separation between the atom and the orientation of the atomic dipole moments with 
respect to the interatomic axis. The analytic expression for reads [18| H9] 

r fi /-^ \2i sin (fcria) 



ri2 = I [1 - (A ■ fuf] 



kru 
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cos (/cri2) sin (kvu) 
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where (x is the unit vector along the dipole moments of the atoms, which we have 
assumed to be parallel {ft = fii = [12), ru is the unit vector in the direction of ri2, 
k = Uq/c, and ri2 is the distance between the atoms. The parameter ri2 is a damped 
oscillatory function of ri2, that varies from ri2 = F for ri2 — )■ to ri2 = for — )■ 00. 
It can have positive as well as negative values. 

We note here that the transient spectrum ([1]) includes effects of a finite bandwidth of 
the detector. In experimental practice, this implies that our calculations are appropriate 
to observations involving a spectrometer (for example, a Fabry-Perot interferometer) 
that analyses the frequency spectrum of the emitted field. 

The two-time correlation function, appearing in Eq. ([T]), is found from the solution 
of the master equation for the density operator of the system and quantum regression 
theorem. We assume that the system we consider is composed of two identical two-level 
atoms (qubits) interacting with a broadband reservoir whose modes are in a squeezed 
vacuum state. We study the behavior of the system in terms of the density operator 
of the atoms, which in the Born-Markov and the rotating-wave approximations satisfies 
the following master equation 

P= j^[Hcs,p]+ Cspp, (3) 

where we have explicitly distinguished two evolution terms, a coherent evolution under 
the non-Hermitian Hamiltonian 

i=l \ jjti 



(2) 
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1 2 

- -th J2 r., {(1 + iV) S^Sr + NSrSi 

- {MS+S^e-^''^'' + M*SiSre^''^^')} , (4) 
and an incoherent evolution that is solely due to spontaneous emission events 



£spP = [(1 + N) SrpS; + NStpS-] 

2 

- [MS^pS+e-^'"^^' + M*5-p5-e2*^=*) . (5) 

Here M = |M| exp(i0s) is the magnitude of two-photon correlations within the squeezed 
vacuum field, (ps is its phase, N describes the number of photons in the field, uq is the 
atomic transition frequency, and Us is the carrier frequency of the squeezed field that we 
shall assume to be equal to the atomic transition frequency. The parameter Q12 = ^21 
stands for the dipole-dipole interaction strength between the atoms [Ml |19] . 

The first line in the effective Hamiltonian describes the coherent evolution 
between the atoms due to dipole-dipole coupling, the second line describes incoherent 
damping of the atoms due to the presence of the broadband squeezed vacuum, and the 
third describes the coherent excitation of the atoms by the squeezed field. The term 
£spP appears as an another source of incoherent damping. 

It is convenient to study the atomic dynamics in terms of the collective states of 
the system, which are the eigenstates of the Hamiltonian 

2 2 

= hY^oS: + ^^JSlSr. (6) 

i=l i^j=l 

Here the dipole-dipole interaction can be viewed as a mixer of the atomic bare states. It 
produces collective states that are found by the diagonalization of the Hamiltonian ([6]): 

|e) = |ei)|e2), \g) = \gi)\g2), 

\s) = -^(|ei)|(72) + |(7i)|e2)), 

\a) = ^{\e,)\g,)-\g,)\e2)), (7) 

where |ej) and \gi) represent the excited and ground state of the zth atom, respectively. 
In the basis of the collective states, the system is composed of four non-degenerate 
energy states, shown in Fig. [H the ground state \g) of energy Eg = 0, two intermediate 
states |s) and \a) of energies Eg^a = h{uo ± 1^12), and the upper state |e) of energy 
Ef, = 2hwQ. The existence of the maximally entangled states, the symmetric \s) and 
anti-symmetric \a) states, refiects the presence of the coherent dipole-dipole interaction 
between the atoms. 

Since atoms interact with each other, the collective state picture provides a good 
approach for studying the initial phase problem. In what follows, we use the collective 
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Figure 1. The energy- level diagram of a two-atom system with the allowed 
spontaneous transitions occurring with rates T + T12 and F — ri2. 



states as the basis for the density operator of the system, and find that the two-time 
correlation function appearing in Eq. ([T]) can be written in terms of the density matrix 
elements as 

2 

J2 r.ASnt)S-{t + r)) = (r + r^^) [Pes{t + r) + p.,,{t + r)] 

+ (r - Tu) [Pagit + r) - p,ait + t)] , (8) 

with the initial r = values determined by the one-time density matrix elements 

Pes(t + T)\r=0 Pee{t) , Psg(t + T)\r=0 Pss(t), 
Pea{t + T)|r=0 - Pee(t) , Pag(t + T)\r=0 -> Paa{t), 

Pgsit + T)\r=0 Pge(t), Pga(t + r)|^=o -Pge(t)- (9) 

One can see that the correlation function appears as a sum of the coherences which 
cascade down the atomic levels. Note that there are no correlations in the transitions 
between and inside the cascade down channels |e) — )■ |s) — t- \g) and |e) — t- \a) — )■ \g). This 
property allows as to analyze spectral properties of each of the transitions separately. 

To introduce initial-phase spectroscopy, we first consider a simplified two-atom 
system, namely the Dicke model, in which the transition channel |e) — )■ \a) — )■ \g) 
does not participate in the dynamics of the system. As we shall see, this will allow us 
to describe and characterise initial-phase spectroscopy using analytic methods. What 
do we mean by "initial-phase" spectroscopy? A squeezed vacuum itself depends on 
the phase 0^, but no phase dependent effects are observed until a reference phase is 
introduced. Usually, the reference phase is provided by a coherent laser field. 

An another way to introduce the reference phase is to prepare the system in an 
entangled state and study the dependence of the evolution on the relative phase A0 
between the initial entangled state and the squeezed vacuum. This is what we call 
the "initial-phase spectroscopy" that could allow us to study the dependence of the 
decoherence process on the phase and to identify entangled states that could be 
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prevented from decoherence. We look for signatures of slowed down or accelerated 
decoherence in the transient spectrum of the emitted fluorescence field. 

Usually, decoherence processes in a given system are studied by analyzing widths 
of the spectral lines. They are determined by asymptotic, exponentially decaying 
correlation functions whose Fourier transforms give spectral (Lorentzian) lines. A 
narrowing of the spectral lines below the natural linewidth is interpreted as a 
stabilization of the quantum fluctuations and hence the reduction of decoherence. 
However, the recent studies on entanglement have revealed that some entangled states 
can decay is a non-exponential way leading to the phenomenon of entanglement sudden 
death [39lllll 1121 HH] . This kind of decay is not reflected in the spectral linewidths. As we 
shall see, the non-exponential decay can be reflected in magnitudes of the spectral lines. 



3. Evolution of the density matrix elements 

The spectrum is determined by the coherence terms that govern spontaneous and 
stimulated transitions in the system. The coherences are solutions of the equations 
of motion that are found from the master equation ([T]). In the rotating frame oscillating 
with the frequency of the squeezed field Ug, the equations of motion for the coherences 
responsible for transitions in the channel involving the symmetric state are in the form 

Pes = -[l^ [^^2 + a) + 1] - ^fiiaj Pes 

+ ^T{n - 1)(1 + a)psg + MT [apgs - (1 + a) , 
Kg = -{l^ [^(2 + a) - 1] + ^fiiaj Psg 

+ (n + 1) (1 + a) Pes + MT [apse - (1 + a) Pgs] , (10) 

and the equations of motion for the coherences responsible for transitions in the channel 
involving the antisymmetric state are 

'Pea = -[l^ [^(2 - a) + 1] + iQu^Pea 

- ^r(n - 1)(1 - a)pag + Mr [apga - (1 - a) Pae] , 

Pag = ~ [^(2 - a) - 1] - ^^12 1 Pag 

+ ^T{n + 1) (1 - a) Pea + MT [apae - (1 - a) Pga] , (11) 

where n = {2N + 1) and a = T12/T is a dimensionless collective damping parameter. 

The coherences oscillate with frequencies shifted from cuq by the amount Q12, the 
dipole-dipole interaction strength, and the squeezed field couples the coherences inside 
the channels. In addition, the squeezed correlation M couples the coherences to their 
conjugates. It is interesting to note that there is no coupling between the symmetric 
and antisymmetric channels. 
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The solution of Eqs. ( TTUj) and ( TTT]) depends, of course, on the initial state of the 
system. We shall assume in all our considerations that the initial state of the atomic 
system is the maximally entangled Bell state 

mt = 0)) = ^{\g)+e^^-\e)), (12) 

which corresponds to an excitation of the two-atom system into a coherent superposition 
of the states \g) and |e). Here, 0b is the (fixed) phase of the state that serves as the 
reference phase for the squeezed field. The application of the phase dependent initial 
state allows the phase control of the fluorescence spectrum and the evolution of the 
initial entanglement. 

The initial non-zero density matrix elements of the system prepared in the entangled 
state (fT2ll . before the application of the squeezed field, are 

Pee(O) = Pgg{0) = ^, 

Pe,(0) = p^,(0) = ie-^-, (13) 

and all the remaining density matrix elements are zero. 

Our task is to trace the time evolution of the spectrum from the initial state 0121) 
to the steady-state distribution. The evaluation of the transient spectrum requires the 
knowledge of the time evolution of the populations of the collective states and the two- 
photon coherence between them. These are determined from the following equations of 
motion 

Pee = -in+l)Tpee + 2\M\Ti2pu 

+ i(n - l)r [(1 + a) + (1 - a) Paa] , 
Pss = ^ (1 + «) r [(ti - 1) - {3n - l)pss 

- {n~l)paa + 2pee-2\M\pu], 

Paa = - (1 - a) r [(n - l)-(3n - l)paa 

- {n- l)pss + 2pee + 2|M|pJ , 
Pu = Ti2\M\ - nTpu 

- \M\T [(1 + 2a) Pss - (1 - 2a) p J , (14) 

where p„ = {peg exp(— 10..) + pg^ exp(i0s))/2 is the real part of the two-photon coherence. 

Note that the initial value of Pu{0) = cos(A0)/2 depends on the relative phase 
A0 = (pB — (ps between the initial entangled state and the input squeezed field. This 
is the quantity that brings the phase dependent features into the time evolution of the 
fluorescence field and its spectral distribution. It is also interesting to note that the 
evolution of the populations is insensitive to the dipole-dipole coupling strength Q12. 
This is due to the broadband nature of the applied squeezed field, which makes the 
squeezing parameters and M independent of frequency. 
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4. Dicke Model 

We first concentrate on a simplified model of the two-atom system, the Dicke model, 
where the evolution of the system through the antisymmetric channel is ignored. 
Physically, it corresponds to a system of two atoms confined to a region which has 
linear dimensions much smaller than the resonant wavelength, ri2 ^ Aq = 2ttc/ujo, so 
that both atoms experience the same amplitude and phase of the field. In other words, 
the atoms act collectively with a common dipole phase. Mathematically, the model 
is described by the same set of equations of motion as Eqs. ffTOj) . ffTTl) and ([H]), but 
with replaced by T. Effectively, the model behaves as a three-level cascade system 
with upper state |e), intermediate state \s) and ground state \g). The transitions in the 
cascade are both damped with the rate 2T. 



s{cjj.)/r 

0.50 r 

0.25 - 








Tt 



Figure 2. (a) The transient spectrum S{uj,t) as a function of the normahzed 
frequency uj — (uj — i-uq)/T and time Tt for the Dicke model with — 2T, A(f> = 0, 
N = 0.5, = 20r, and = y/N{N+l). The frame (b) illustrates the time 
evolution of the populations pee (thick line) , and pss (thick dashed line) together with 
the two-photon coherence p„ (thin line). The atoms were initially in the Bell state (|12p . 

In what follows, we calculate the transient fluorescence spectrum by solving the 
equations of motion for the density matrix elements with = F, from which we then 
evaluate the spectrum ([T]) in the the limit of well separated spectral lines, ^ TN. 

Figure |2t^a) shows the transient spectrum calculated for the relative phase A(f) = 0. 
The system was initially prepared in the Bell state ([T2|) . We see that the amplitude of 
the peak a.t u = uq — is practically zero for all times, whereas the amplitude of the 
other peak is large. The amplitude of the peak at u = Uq + fli2 increases in time with 
the rate F, it passes through a maximum and then slowly decays towards a stationary 
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value. A close look at the transient evolution of the amplitude reveals the existence of 
a minimum (saddle point) that appears at a time near Ft = 1. The existence of the 
saddle point can be interpreted as a result of switching the atomic dynamics from that 
governed by the initial correlation to that governed by the correlation transferred to the 
atoms from the squeezed field. 

In order to interpret the numerical results, we calculate analytically the transient 
spectrum in the limit of Q12 ^ TN corresponding to the case of well separated spectral 
components. The spectrum depends on the bandwidth of the detector, for simplicity 
and without loss of generality we will assume that F^^ — )■ 00. The limit of F^^ — )■ cxd 
corresponds to a very broad detector. In this case, the spectrum ([T]) reduces to a finite 
Fourier transformation of the two-time atomic correlation 

2 

S{u,t) = Re V Tij / dTe-"^^{Slit)Sr{t + r)). (15) 

The spectrum might be interpreted as the rate of production of photons of frequency oj 
at time t. 

Using Eqs. (fTO!) and (fTTi) . and neglecting all terms of order ri^2 smaller, we 
obtain the following expression for the spectrum 

2{n + l)pUt)-2\M\pu{t) , (n - l)pee(t) + |M|p4t) 



Siu,t) = 2F 



+ Pss{t) 



2{n-l) ^ n+1 



(16) 



where \M\ = \/n^ — 1, a) = (w — a;o)/F and ^\2 = ^12/F. The spectrum has a simple 
structure. In general, at a given time, the spectrum is a sum of four Lorentzian lines, 
with two overlapping lines at each of the transitions. The amplitude of the lines at 
the transition frequency u = Uq — Q12 is proportional to the population Pss{t) of the 
symmetric state, whereas the amplitude of the lines at u = Uq + fli2 is determined by 
the population peeit) of the upper state. However, most important from the point of 
view of this paper is that there is an additional contribution to the amplitude of the 
line at CO = uq + Qu from the coherence Peg(t). Evidently, the two-photon coherence 
and in particular its dependence on the relative phase can directly affect spontaneous 
transitions originating from the state |e). 

It is clear from ( IT6l) that the absence of the peak at u = uo — Q12 results from the 
absence of the population in the symmetric state. This is also seen from Fig. Mjo), where 
we plot the time evolution of the populations Peeif), Pssif) and the coherence puit). The 
absence of the population in \s) is readily understood if we recall that there are two 
pathways for excitation of state |e). First, there is a stepwise excitation in which the 
system goes from \g) to \s) and then from \s) to \e). Second, there is a direct transition 
from \g) to \e) involving a simultaneous absorption of two photons. The two-photon 
correlations, the initial and that contained in the squeezed field, cause simultaneous 
absorption of two photons which results in the direct transfer of population from the 
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ground state \g) to the upper state |e) without populating the intermediate state \s). 
Thus, the buildup of the peak at w = — ^^12 is suppressed. 



0.50 r 

0.25 - 




rt 

Figure 3. Same as in Fig. [21 but for the relative pliase — tt. 

The dependence on time of the spectral features changes qualitatively when A0 is 
varied. A change of A(j) from to tt leads to a hole burning in the higher frequency 
peak. This is shown in Fig. [3l where we plot the spectrum for the same parameters 
as in Fig. |2] but A0 = vr. We observe a rapid collapse of the amplitude of the peak 
at a; = Wo + ^^12, which then evolves through a hole burning to a narrow peak as time 
varies from t = to t = 2r. To explain the origin of the hole burning, we illustrate in 
frame (b) the time evolution of the populations Peeif), Pssif) and the coherence Puif)- 
Evidently, at time when the hole is burned in the spectrum, the population pee{t) is 
almost zero and Pu(t) is negative. 

The hole burning, as illustrated in Fig. [3]^a) and its persistence for only small 
period of time around Ft = 1 can be understood as the result of the destructive 
interference between the initial and squeezing correlations that oscillate out-of phase 
when A(f) = n |55|- An alternative explanation is as follows. The spectral line at 
u = uq + Q12, as seen from Eq. ( 1161) is composed of two Lorentzians of different widths 
and weights. Each of the weights can be positive or negative depending on the sign of the 
coherence Pu{t)- For A0 = 0, the coherence pu{t) is positive for all time and much larger 
than the population peeif)- In this case, the weight of the broader spectral component 
with the width, 2nT is negative while the weight of the narrower spectral component 
with the width, nF, is positive. MoUow |56] and Rice and Carmichael [57] pointed 
out that the negative weight of the broader Lorentzian can give rise to a spectrum 
with subnatural linewidth. Thus, the subtraction of the broader Lorentzian from the 
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narrower one results in a narrowing of the spectral line and the reduction of the intensity 
of the line. The narrowing of the spectral line is not well visible in Fig. 2, but the 
reduction of the amplitude, a saddle point near F/f: = 1 is well visible in the figure. The 
situation changes dramatically when one changes the relative phase to A0 = vr. In this 
case, Pu{t) is negative for short times so that the narrower Lorentzian can now have 
a negative weight while the broader Lorentzian has a positive weight. Thus, for times 
close to the value which minimizes the weight of the narrower Lorentzian, we obtain a 
hole burned in the spectrum [57] . 

The enhancement of the spectral peak at u = Uq — Q12 reminds us the phenomenon 
of the Dicke superradiant pulse formulation [581 159] . Over a short time period, the 
population is transferred to the symmetric (superradiant) state and is then emitted with 
the enhanced, superradiant rate F + F12, which is the damping rate of the population 
of the symmetric state. 

The hole burning can be connected with the phenomenon of entanglement sudden 
death [391 HH US], HH] . To show this, we compare the time evolution of the spectrum with 
the evolution of an initial entanglement between the atoms. We quantify entanglement 
by calculating the concurrence, which can be evaluated explicitly from the density 
operator of the system as [60] 



where p* denotes the complex conjugate of p and ay is the Pauli matrix. Concurrence 
varies from C = for separable qubits to C = 1 for maximally entangled qubits, and the 
intermediate cases < C < 1 characterize a partly entangled qubits. 

The concurrence is specified by the density matrix of a given system and thus is 
evaluated from the knowledge of the density matrix elements. For the system considered 
here of two atoms interacting with a squeezed vacuum field, the density matrix written 
in the basis of the collective states ([7]) has the X-state form 




(17) 



where the quantities Aj are the eigenvalues in decreasing order of the matrix 

R = p{ay® ay) p* {ay ® ay) , 



(18) 



/ Pgg Pge ^ 

Pss 

Paa 

\ Peg Pee J 



(19) 



This simple form of the density matrix leads to the concurrence of the form 
C{t) =max{0, Ci{t), C2(t)}, 



(20) 



where 



Ci{t) 



2|Pe9(t)| - [Pssit) + Paa{t)] , 



(21) 



and 



C2{t) 



Pss{t) - Paa{t) \ - 2^ Pgg pee- 



(22) 
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It is clear that the concurrence C{t) can always be regarded as being made up of the 
sum of nonnegative contributions Ci{t) and C2{t) that are associated with two different 
sources of entanglement in a two atom system. The contribution Ci{t) provides a 
measure of an entanglement produced by the two-photon coherence Peg, whereas C2(t) 
provides a measure of an entanglement produced by the populations of the maximally 
entangled states \s) and \a). 

It is easily verified, that with the initial Bell state f|T2|) . C2{t) is always negative, so 
that a possible entanglement in the system is solely determined by Ci{t). 



0.5 








Ft 

Figure 4. Time evolution of the concurrence for the Dicke model coupled to a 
broadband squeezed vacuum field with N = 0.5, ~ y/ N{N + 1), and for different 
relative phases: A(f) — tt (solid line), A0 — (dashed line). 

Figure H] shows the time evolution of the concurrence for the same choice of 
parameters as in Figs. 2 and 3. We see a strong variation in the time evolution of 
the concurrence with the relative phase A0. For the case of A(f) = 0, the atoms are 
entangled for all times with the initial entanglement decaying slowly towards a non-zero 
steady state value. The evolution of the initial entanglement is completely different 
when the phase is changed to A0 = tt. In this case, the initial entanglement decays in 
a non-exponential fashion and disappears at a very short time of the evolution. Beyond 
that time, the atoms remain disentangled for a finite period of time, and then the 
entanglement reappears again [iQl HI], and evolves asymptotically to a steady-state 
value that is independent of the phase A0. A close look at the spectrum and the 
concurrence reveals that the entanglement disappears for the period of time when the 
hole occurs in the system. Thus, the hole burning in the spectrum can be regarded as 
a signature of the phenomenon of sudden death of entanglement. Hence, observation of 
the hole burning in the spectrum would be a striking confirmation of the entanglement 
sudden death. 
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5. Non-Dicke Model 

We now turn to the problem of phase effects in an extended atomic system. We relax 
the small sample approximation and consider dynamics of an extended atomic system, 
the so-called non-Dicke model, where we explicitly take into account a non-zero spatial 
separation between the atoms. The separation will act as a phase difference between 
the oscillating dipoles. As we shall see, this will introduce an additional phase to the 
dynamics of the system. In other words, the transient spectrum depends strongly on 
whether > or < 0. 
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Figure 5. (a) The transient spectrum S{(jj,t) as a function of tlie normalized 
frequency cj = — loq) /T and time Tt for the non-Dicke model with = 2r, Ac/) = 0, 
N = 0.5, \M\ = y/N{N + l), = 20r, and Tia = 0.5r. The frame (b) illustrates 
the time evolution of the populations pee (thick line), pss (thick dashed line), and paa 
(dashed-dotted line) together with the two-photon coherence p„ (thin line). The atoms 
were initially in the Bell state (fT2|). 



For phase independent initial states and/or absence of the squeezed field, the 
spectrum of the field emitted from a two-atom system is composed of two peaks, a 
broad peak of the width F -|- F12 due to spontaneous emission from the symmetric state 
and a narrow peak of the width F — F12 arising from spontaneous emission from the 
antisymmetric state |1HI EI]- In this case, the collective damping imposes quantitative 
changes to the spectrum by modifying the widths of the spectral lines, but does not 
lead to spectra which are qualitatively different from those that occur in the absence of 
the collective damping. Here, we show that the inclusion of an initial phase can lead to 
transient spectra which are quite distinct to those seen in its absence. 

The transient spectrum can be found analytically by solving the equations of 
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motion fllUp and fill I) . However, the analytical form of the spectrum is very complicated 
and difficult to interpret. Therefore, we restrict ourselves here to the graphical 
illustration of the spectra. 

Figure [5](a) shows the transient spectrum for A0 = and > 0. The spectrum 
is composed of two asymmetric peaks, whose initial amplitudes decay exponentially in 
time and attain their stationary values after a short time t ~ F"^. The origin of the 
peaks and their time evolution is easily explained from the analyses of the populations 
of the collective states, shown in Fig. ©(b). Since /3ee(0) = 1/2 and Pss(O) = Paa(O) = 0, 
it is clear that the peaks arise from the |e) — > \s) and |e) — )■ \a) transitions. However, 
the spontaneous transitions |e) — )■ |s) do not lead to population of the state \s) that 
remains almost unpopulated for all times. A transfer of the initial population occurs 
only from |e) to \a). Moreover, the transfer of the population to the state \a) occurs on 
the time scale much shorter than t = (F — Fi2)~^, the time scale required for appreciable 
changes of the population of \a) to occur, as predicted by the rate equations (IT^ . This 
indicates that the state \a) is populated by stimulated rather than spontaneous emission 
process. 



s{oj,t)/r (a) 

0.30 r 




i.O 



Figure 6. Same as in Fig. [5] but for the relative phase — n. 

The situation changes when we switch the relative phase from to vr. In this case, 
a hole burning occurs in the higher frequency peak, and simultaneously the amplitude 
of the lower frequency peak becomes enhanced. This is demonstrated in Fig. [6l where 
we plot the spectrum for the same parameters as in Fig. [5], but with A0 = vr. The 
origin of this effect is the same as discussed above for the Dicke model. The explanation 
again follows from the observation that the destructive interference between the initial 
two-photon coherence and that induced by the squeezed field leads to spontaneous 
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processes which dominate over the stimulated processes. As the result, the symmetric 
state becomes significantly populated on a time scale t = (F + ri2)~^, which is the time 
scale for the population of the symmetric state, predicted by the rate equations. This 
shows clearly that the phase A0 changes the transient populations of the collective states 
of the system. This fact is illustrated in Fig. Mjo), where we show the time evolution 
of the populations and the two-photon coherence. Clearly, at the time when the hole 
burning occurs, the symmetric state is significantly populated and the states \a) and |e) 
are essentially completely depopulated. 



S{Cj,t)/r (a) 
0.30 r 




n 

Figure 7. Same as in Fig. [5]but for a negative collective damping rate, = — O.SF. 

It is interesting that hole burning can occur not only by changing the relative 
phase A0 , but also by changing the sign of ri2. A change of the sign of from 

> to ri2 < is achieved by varying the distance between the atoms. This is 
shown in Fig. [TJ where we plot the spectrum for the same parameters as in Fig. |5l but 
for Fi2 = — O.SF. We see that the sign of the collective damping parameter, which 
depends on the phase difference between the oscillating atomic dipoles, has a strong 
impact on the shape of the spectrum. A change F12 — > — ri2 leads to hole burning in 
the lower frequency peak. Thus, the effect of negative F12 is seen to be the same as the 
relative phase A0 = vr. This is a surprising result as one could expect that a change of 
the sign of ri2 would result in a spectrum similar to that as shown in Fig. |5l but with 
the peaks exchanged in position. The changes in the spectral distribution indicate that 
the phase of the squeezed field is locked to the phase difference between the oscillating 
dipoles. We may conclude that the collective damping acts like an artificial tweezer that 
rotates the phase of the squeezed field. To say this another way: The squeezed field 
keeps memory of the relative phase of the atomic dipole moments. 
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Figure 8. Same as in Fig. [5] but for Acf) = tt and T12 = — O.SF. 

It is not difficult to see that may act as a tweezer of the squeezing phase. This 
fact is graphically shown in Fig. [HI where we plot the transient spectrum for the same 
parameters as in Fig. [TJ but for = tt. When A0 is switched from zero to tt, the 
spectrum reverts to a profile similar to that shown in Fig. [5]^a) that the spectrum is 
composed of two asymmetric peaks. 

The role of as a phase shifter is also seen from the solution for the two-photon 
coherence peg- It is not difficult to show from Eq. ([T3j) that the long time [t — t- 00) 
solution for the two-photon coherence Peg is of the form 



Peg = 





a||M|e^('^»+*'") 


n2[n4 + 4|M|2(|a 


2 _ ^2)] 



(23) 



where is a phase angle with only two discrete values of $a = for > and = tt 
for < 0. It is clear that a change ri2 —T12 leads to a shift of the squeezing phase 
by TT. It is also interesting to note that the magnitude of the coherence is proportional to 
ri2|M|, which indicates that the collective damping ri2 can be regarded as an efficiency 
of the transfer of the two-photon correlations \M\ from the squeezed field to the atoms. 

Finally, we compare the transient spectrum with the transient properties of the 
concurrence. Figure [3 shows the concurrence C{t) as a function of time for the same 
parameters as in Figs. [5] and [71 The decay time of the initial entanglement depends 
on A(f). We observe that similar to the Dicke model, the initial entanglement collapses, 
i.e. disappears over a finite time but, what is surprising, it never revives. The lack of 
recurrence of entanglement in the non-Dicke system is due to the fact that the presence of 
the antisymmetric state increases the threshold for entanglement in Ci{t) from pss{t) to 
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Figure 9. Time evolution of the concurrence for tlie non-Dicke model coupled to a 
broadband squeezed vacuum field with N = 0.5, \M\ — y/ N{N + 1), Acj) = 0, and 
different signs of ri2: = — 0.5r (solid line), — O.SF (dashed line). 

Pss(t) + Paa(t) , as seen from Eq. fl2T]) . and decreases the positive term in C2(t) from Pssit) 
to \pss(t) — Paait)\, as Seen from Eq. (1221) . The enhancement of the threshold in Ci(t), 
it turns out, is sufficient to wipe out any entanglement. Physically, the lack of the 
entanglement revival is the result of an imperfect transfer of the two-photon correlations 
from the squeezed field to the atoms. It is easy to see from Eqs. iHM and ( |23l) that the 
two-photon coherence peg is proportional to a = ri2/r that is smaller than one (a < 1) 
for the non-Dicke model. 

6. Conclusion 

We have examined the transient spectrum of the radiation field emitted by a system of 
two atoms that interact with a broadband squeezed vacuum field and initially prepared 
in a maximally entangled Bell state. The initial entangled state provides a reference 
phase for the phase dependent squeezed vacuum field. We have considered the Dicke and 
non-Dicke models. The non-Dicke model explicitly includes a phase difference between 
the atomic dipole moments. Our calculations have demonstrated the qualitatively new 
behavior observable in the transient spectrum of a collective system that evolves from an 
initial entangled state. It has been found that the effect of the relative phase suppresses 
one of the two spectral peaks normally expected in transient spectra and enhances the 
other. The suppression is seen as a hole burning in the spectral line. The physical origin 
of the hole burning has been explained in terms of the destructive interference between 
the initial and the squeezed field induced two-photon coherence. The relation between 
the phases determines whether the population is transferred to the two-excitation energy 
state or to a single-excitation collective state. We have also made a comparison of the 
transient behavior of the spectrum with the time evolution of the initial entanglement 
and found that the hole burning can be interpreted as a manifestation of the phenomenon 
of entanglement sudden death. In addition, we have found that in the case of the non- 
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Dicke model, the collective damping rate that can be positive or negative, introduces a 
phase shift of the phase of the squeezed field. In all cases we have considered, we have 
found that a change of the sign of the collective damping results in the hole burning in 
one of the spectral lines, and the suppression of a given peak determines the state that 
is selectively populated. 
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